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Abstract 

When an A-photon state enters a lossless symmetric beamsplitter from one input port, 
the photon distribution for the two output ports has the form of Bernouli Binormial, with 
highest probability at equal partition (N / 2 at one outport and N / 2 at the other). However, 
injection of a single photon state at the other input port can dramatically change the photon 
distribution at the outputs, resulting in zero probability at equal partition. Such a strong 
deviation from classical particle theory stems from quantum probability amplitude cancella- 
tion. The effect persists even if the A-photon state is replaced by an arbitrary state of light. 
A special case is the coherent state which corresponds to homodyne detection of a single 
photon state and can lead to the measurement of the wave function of a single photon state. 


1 Introduction 

Interference effect of light has played an important role in the conceptual development of quantum 
theory. Richard Feynmann once wrote 1 that the Young’s double slit experiment “has in it the heart 
of quantum mechanics”. But the phenomena of interference do not simply stop at Young’s double 
slit experiment and its variations. Much richer phenomena occur in higher-order interference 2-6 
when there are more than one particle involved in the process. For example, Greenberger et al. 7 
recently proposed new demonstration of locality violation by quantum theory with superposition 
state of three or more particles. 

In the meantime, along a quite different line, Ou and Mandel 8 have investigated a startling 
quantum interference effect where a strong field interferes with a considerably weak field. It was 
shown 3 and demonstrated 8 that for certain nonclassical fields, the interference fringe visibility does 
not change even though the ratio of the intensities of the two interfering fields is much greater 
than 1, in conflict with the intuitive picture from classical wave theory for interference. In this 
case, the seemingly insignaficant weak field plays an essential role for the interference effect even 
though its intensity is negligibly small. Therefore, the presence of the weak field can dramatically 
change the outcome of the result. 

In this paper, we will present another example of how existence of a weak field can make 
a significant difference. It deals with N + 1 photons with N being a positive integer. We will 
consider a situation when an N-photon state interferes with a single photon state with the help 
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of a symmetric lossless beamsplitter (see Fig.l). A special case of N = 1 has been experimentally 
investigated as an example of fourth-order interference. 9 However, quite different from the two- 
photon coincidence measurement technique used in fourth-order interference, we will exam photon 
probability distribution at two output ports of the beamsplitter. Although no interference pattern 
exists, the phenomenon discussed here attributes to quantum interference of multi-particle ( N + 1 
particles). We will also extend the discussion to an arbitrary state input in replacement of the 
N- photon state. 



FIG. 1. Layout for the interference between N -photon state and a single photon 
state via a beamsplitter. 


2 Photon Probability Distribution for a Symmetric Loss- 
less Beamsplitter 

It is well-known that when a number of particles, say N, enter a 50:50 lossless beamsplitter from 
one input port, the particles are randomly sent to the two output ports with equal probability, 
resulting in the simple Bernoulli binormial distribution as 

N\ 

Pq(Ni,N 2 ) = 2 n n \n 2 \ C 1 ) 

Ni is the number of particles exiting from output port 1 while N 2 is for port 2. In the case 
of photon, the above result suggests that each photon acts independently as a classical particle. 
The wave behavior of light does not show up here because of the absence of superposition. This 
distribution has its maximum when = N 2 = N/2 (equal partition). So it is most likely to find 
equal number of photons on each side of the beamsplitter. For large N and \N\ — A^| << Af, 
Eq.(l) becomes 

Po{N u N 2 ) = -£= 6 Nl+NliN (2) 

V2Nic 

which is a Gaussian. The extra factor of 2 is because Po(Ni, N 2 ) = 0 for every other value of 
- N 2 . 

Next, we let a single photon state enter the input port 2 of the beamsplitter. We will look 
for the probability P\(N\, N 2 ) that N\ photons exit at output port 1 while the other N 2 photons 
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at port 2 with N\ + N 2 = N + l. Let us for a brief moment consider the outcome from classical 
particle theory. As a classical particle, the input single photon will have 50% of probability going 
out at either ports. Because the single photon is independent of the other N photons, we simply 
add the probabilities to obtain the final result: 


P?{Ni,N 2 ) 


\N\ \N\ 

2 2 N (Ni - 1)\N 2 ] + 2 2 n (N 2 - 1)!A\! 


(N + 1)! 
2 N+l N l ]N 2 \' 


(3) 


which is in the exactly same form as that in Eq.(l). Therefore the existence of the single photon 
at the other port does not influence the photon probability distribution at all. The single photon 
from port 1 acts as if it were part of the N photons from the port 1. This is because classical 
particles are independent of each other and it doesn’t matter which port it enters. 

On the other hand, the outcome is totally different if we treat the photons as quantum particles. 
We cannot simply add the probabilities. The principle of quantum mechanics requires that the 
probability amplitudes be added. For simplicity, let us first consider the case when N is an odd 
integer and jVj = N 2 = (N + l)/2. The probability amplitude has two contributions as shown in 
Fig. 2: (a) the single photon input at port 2 goes directly to output port 2 while JVj — 1 = (N — 1)/2 
of the N photons input at port 1 are reflected and go to output port 2 and N 2 = ( N + 1 )/2 photons 
to port 1, or (b) the single photon is reflected and goes to output port 1 while N 2 — 1 = (N — l)/2 
photons go to output port 1 and N 2 = (N + l)/2 photons are reflected to port 2. From Eq.(l), 
we find that these two possibilities have equal probability thus their probability amplitudes have 
equal absolute value. For their phases, however, because there is a 7r/2 phase shift for the reflected 
field and no phase shift for the transmitted one at a symmetric beam splitter, the total phase shift 
for the N 1 photons at the output ports will be different for the two possibilities. Referring 
to Fig.l, we find that the total phase shift for the first possibility mentioned above is tp a = 
(N 2 - 1)tt/2 = (N - 1)7t/ 4 while for the second possibility, (p b = tt/2 + N 2 n/2 = (N + 3)7t/ 4. The 
phase difference between the two possibilities is thus <p b — Va = Therefore, the two probability 
amplitudes will cancel each other, resulting zero probability for N x = N 2 = (N + l)/2. This 
result is completely different from that of a classical particle theory in Eq.(3). As seen above, the 
probability cancellation at Ni = N 2 results from the quantum interference of N + 1 particles. 



N+l 

2 



_ N+l 
2 


FIG. 2. Two contributions to the output photon distribution. 
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For the other cases when N\ 7 ^ N 2 , we cannot use the simple argument as above. But we may 
derive the output state along the line of Ref .10 and find the probability distribution P\(N U N 2 ). 
Or we can use the formula 

P,(NuN 2 ) = (: ( 4 ) 


A\ — (®i + * o, 2 )j \/2, A 2 = (a 2 + i o,\ )/y / 2 

are the annihilation operators for the output modes for a symmetric lossless beamsplitter. The in- 
put modes represented by a x ,a 2 are in the state of |$) = |/V) 1 | 1 ) 2 . After some lengthy calculation, 
we have 

AH 

P l (N l ,N 2 ) = 2 /v+i yVj!/v 2 ! ~ N2 ^ 6n ' + /v 2 - /v +i- ( 5 ) 

The above expression can also be derived from the general formula given by Campos, Saleh and 
Teich in Ref.il for arbitrary numbers {n x ,n 2 } of input photons at the two input ports, with the 
setting of r = 1 / 2 , n x = N,n 2 = 1 . When N,N X ,N 2 » 1 , Eq.(5) can be approximated by 

W,^)« ( *' l i, iVl)1 A, e- ,N, - A ' ,|1/2 "W«„N + .. (6) 


2 

\/ 2 ~Nn 


Notice that when N is an odd integer, P\(N X ,N 2 ) — 0 for N x — N 2 = (N + 1 )/ 2 , exactly as 
predicted from the simple argument of probability superposition given in the previous paragraph. 
When N is an even integer, P x (N/2 + \,N/2) = N\/2 n+1 (N/2 + 1 )!( A^/2)! ^ 0, but because 
P x (N/2 + 1 , N/2)/ P 0 (N/2 + l,A r /2) = l/(N + 1 ) << 1 for N >> 1 , or the probability with 
a single photon input is much smaller than that with vacuum state input, the probabilities for 
Ni « N 2 are quite different in the two cases with or without the single photon state at port 2. 
Actually, the whole probability distribution in Eq.(5) is different from the probability distribution 
in Eq.(l), as seen in Fig. 3. The maximum probability for Pi(N x ,N 2 ) occurs at \N\ — N 2 \ Vn 
or N x tt(N± s/N)/2 while for Po^, N 2 ) it occurs at N x = N 2 rs N/2. The existence of a single 
photon dramatically changes the pattern of the output photon distribution. 



N, N, 

FIG. 3. Output photon distribution for N-photon state input at port l with (a) 
vacuum state or (b) single photon state at port 2 (N=19). 
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3 Interference of a Single-Photon State with Arbitrary 
State 


The above quantum probability cancellation effect due to a single photon state is not stricted to 
/V-photon state as input state. Let us consider an arbitrary state of light input at port 1. Its 
state is generally described by the Glauber P-distribution P tn (a). But before going into lengthy 
calculation, we may take a guess about the photon distribution of the output fields by the following 
argument: since the vacuum state and the single photon state are completely incoherent in the 
sense that they have a totally random phase distribution, the output fields due to interference 
of one of these states with any other state will not have any coherence information of the input 
state. Therefore, the output photon distribution of the beamsplitter will lose all the coherence 
information of the input state and will depend simply on the photon statistics P^ n of the input 
state at port L So combining this fact with Eqs.(l,5), we come up with the output photon 
distributions in the form of 


Po(N u N 2 ) 


W + tt.)! pi, 

2 n '+ n i Ni\N 2 \ rNl +Ni 


(7a) 


for vacuum input at port 2 and 


Pi(N u N 2 ) 


(N x + N 2 - 1)! 

2 n ' +n * Ni\N 2 \ 


( N , - /V ,) 1 


(76) 


for single photon state input at port 2. Of course, we may rigorously derive the output photon 
distribution by following the procedure described in Ref. 10 to first find the state of the output 
fields of the beamsplitter in terms of the P-distribution. The photon distribution for the output 
fields can then be calculated through Eq.(4). It can be shown that Eq.(7) is indeed the correct 
form for the output photon distribution. 

By comparing Eqs.(7a) and (7b), we easily find that P\(N\ = N 2 ) = 0 for single photon state 
input at port 2 while 

p ( _ n \ _ y' ( 2 ^)- ^ q 

for vacuum input. Therefore, the existence of the single photon state at port 2 does make a 
difference in the output photon distribution even for arbitrary input state at port 1, and the 
probability for N x = N 2 is exactly equal to zero. The cancellation of the probability for jVi = N 2 
is because of the destructive interference between the N photons and the single photon as we 
discussed above. 

In an actual experiment, however, it is difficult to measure the complete distribution P(N X , N 2 ), 
but the distribution P(N X - N 2 = M) can be measured by balanced homodyne detection. 12,13 From 
Eqs.(7a,b) we find that 


Po{N x -N 2 = M) 


V W - M)! n 
^ 2 2Nl ~ M N X \(N X - M)\ 2N '~ M 


Ni =M 


Pj( N X -N 2 = M) 


OO 

M 2 £ 


N\ =M 


(2M — M- 1)! pin 

2 2N i~ m N X \(N X — M)\ 2N ‘- M -' 


( 8 ) 
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for M > 0. For M < 0, the symmetry between N X ,N 2 in Eq.(7) leads to P(M) = P(-M). 

Next, we will evaluate P 0 {M), P X (M) for some special states. For iV-photon state input with 
N » 1 , we have P* n = 8 n> N, and Eq.(8) gives results similar to Eqs(l,5): 


Po(M) = 


N\ 


2 N (N/2 + M/2)\(N/2 - M/2)! 


V2N 


-m 2 /2N 


7T 


for N >> 1 


Pi(M) = 


M 2 N\ 


(9) 


2 N + l {N/2 + M/2 + \/2)\{N/2 - M/2 + 1/2)! 

2 M 2 

** 1 W - L nsi * I't r\j * . r . 

for N » 1. 


v/2A^ N 


,-M 2 /2N 


For coherent state input, P£ = n"e "/n! with n being the average photon number. Therefore, we 
have 


°° (9 N, + MV -IN-i+M -n 

P °^ = 2 27Vl+M ^V 2 !(iV 1 + M)! (2 yv 2 + M)! = 6 " /yif ^ 

p fM) y* M 2 {2N\ + M — 1)! ^N 2+ M-i e -n M 2 ,_ v 
1 2 2/Vl+M jY 2 !(A f 1 + M)! (27V 2 + M — 1)! — fi 6 

where /m(«) is the Bessel function with purely imaginary argument and has the form of 

/ M (n) = r d<pe~ iM,fi e ncosv « -^1= e*- M2/2fi when n » 1. 

J \/2mr 


Therefore, for large n, 


( 10 ) 


( 11 ) 


Po(M) 


Pi(M) 


V2 


a -M*/2n 


me 


1 M 2 




,-M 2 /2n 


n7r n 


(12) 


Eq.(12) has the same form as Eq.(9) for large N besides the factor of 2 which is explained earlier 
right after Eq.(2). This is not surprising if we consider the fact that when the photon number is 
large, the interference scheme discussed above becomes homodyne detection scheme. Since both 
vacuum state and single photon state have random phase distribution, homodyne detections with 
N-photon state ( N >> 1) and coherent state as local oscillators are equivalent. As a matter of 
fact, the output photon distributions will always have the form of Eq.(12) for any state as local 
oscillator, provided that the average photon number is large and photon number fluctuation is 


much less than average photon number (^/(An 2 ) << n). We can see this point from Eq.(8): when 
y(A n 2 ) << h, Py" has a narrow peak around n and is a fast changing function as compared with 
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other terms in the summation, therefore the contribution to the summation is only from the few 
terms around n, so that we can pull all other terms out of the sum, that is, 


Po(M) 


n\ 


2*(n/2-M)!(n/2 + M) 1 


E 02 


V5 


; M2 / 2n when n >> 1, 


n7r 


(13a) 


and similarly 

1 A/ 2 

P\(M) ~ -7==— e~ M ' 2n when n >> 1. (136) 

\/2n7r n 

We can also understand this result from the fact that any fluctuation in local oscillator is cancelled 
in balanced homodyne detection scheme. 12 

Furthermore, if we set n — ► 00 , we can replace the discrete variable M with a continuous 
one defined by x = M / \fh and the probability distributions in Eqs.(13a,b) lead to probability 
densities of continuous variable x as 


Po(x) 



Pl(x) 



(14) 


which correspond to the square of the absolute value of the wavefunction for the ground state 
and single photon state, respectively. Thus by measuring P(M) in homodyne detection, we can 
deduce the wavefunction of the input state at port 2. This is exactly the technique of optical 
tomography used by Smithey et al. 13 But here we applied it to a single photon state (input at 
port 2) and proved that the outcome does not depend on the state of the local oscillator (input 
field at port 1) as long as the average photon number is large and the fluctuation is not very large 
for the local oscillator (i.e., the condition for the approximation in Eqs.(13a,b)). 

However, there is an exception to the above. It is well-known that for thermal light, we have 


(An 2 ) = n(n + 1) 


so that ^(An 2 ) « n and we cannot use the approximation in Eqs.(13a,b). For thermal light, 
P™ n — h n /(fi + l) n+1 , so from Eq.(8), we have 


Po(M) 


P i(M) 


E 


(2N 2 + M)\ 


n 


2 n 7 +m 


n 2 -- 


0 2 2 ^+M(jV 2 + M)\N 2 \ (n + l) 2 Af 2 +M+i 


M 


n -f 1 


n 


M + 1 M 


2 1 2 

(2iV 2 + M- 1)! 


+ l,M + l;4x 2 ) 


(15) 


n 


2N-2+M-1 


M 2 2 2N > +m (N 2 + M)\N 2 \ (n + 


Mx 


M 


rt 


-H 


M + 1 M 

2 ’ 2 


,M+ l;4x 2 ) 


where x = n/2(n + 1) and F(a, 0,r, z ) is the hypergeometric function. With some re-arrangement, 
we can prove that Eq.(15) have a simpler form as 
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Po(M) 


1 


Pi(M) = 

with 9=1 + 1/n — \/2n + 1/n. For 

Po(M) = 

P X (M) = 

Therefore, The output photon distribution for thermal light input is different from that of coherent 
state input. But the general trend in the change of the shape from P 0 {M) to P\{M) is similar in 
both states (Fig.4). The quantum interference effect due to single photon is the same. 
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FIG. 4 * Probability distribution P 0t i(M) for the balanced homodyne detection of 
vacuum state and single photon state with (a) coherent state or (b) thermal state as 
local oscillator, n = 300. 
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\/2n + 1 

M 


(M > 0) 


(16) 


M 


n 


large n, q M becomes e so that Eq.(15) is changed to 


1 0 -M/y/Z/2 


\j2n + 1 


E e -M /v A7i 

n 


(M > 0) 


(17) 
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It is interesting to note that the weak nonclassical state (single photon state) plays an important 
role in the interference with a strong classical field (coherent state or thermal state) in contrast to 
the case discussed in Ref.3 where the nonclassical interference occurs between a strong nonclassical 
field and a weak classical field. Even though the nonclassical field is weak here, the result is very 
nonclassical in the sense that the probability of detecting equal intensities in the two outputs is 
zero ( Pi(M = 0) = 0). It can be proved that in the similar situation (one field is weak and the 
other is strong), classical wave theory predicts that the probability is largest for equal intensity 
output at the two ports. 

So far we have only discussed the single mode situations. In practice, we always have wide 
spectrum. Since two different sources of light are involved in the interference, the observation of 
the probability cancellation effect requires the overlap of both spatial and temporal mode structure 
of the two fields as well as near unit quantum efficiency of the detectors. 
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